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Neutrinos can play an important role in the evolution of the Universe, modifying some of the cosmo- 
logical observables. In this contribution we summarize the main aspects of cosmological relic neutrinos 
and we describe how the precision of present cosmological data can be used to learn about neutrino 
CN properties, in particular their mass, providing complementary information to beta decay and neutrino- 

less double-beta decay experiments. We show how the analysis of current cosmological observations, 
. ^ such as the anisotropics of the cosmic microwave background or the distribution of large-scale struc- 
^ ture, provides an upper bound on the sum of neutrino masses of order 1 eV or less, with very good 
^ perspectives from future cosmological measurements which are expected to be sensitive to neutrino 
masses well into the sub-eV range. 



0.1 Introduction 

The subject of this contribution is the role of neutrinos in Cosmology, one of the best examples of 
the very close ties that have developed between nuclear physics, particle physics, astrophysics and 
cosmology. Here we focus on the most interesting aspects related to the case of massive (and light) 
relic neutrinos, but many others that were left out can be found in [1, 2j. 

We begin with a description of the properties and evolution of the background of relic neutrinos 
that fills the Universe. Then we review the possible effects of neutrino oscillations on Cosmology. 
The topic of neutrinos and Big Bang Nucleosynthesis is reviewed in a different contribution to this 
Special Issue [3] . The largest part of this contribution is devoted to the impact of massive neutrinos on 
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cosmological observables, that can be used to extract bounds on neutrino masses from present data. 
Finally we discuss the sensitivities on neutrino masses from future cosmological experiments. 

Note that massive neutrinos could also play a role in the generation of the baryon asymmetry of 
the Universe from a previously created lepton asymmetry. In these leptogenesis scenarios, one can also 
obtain quite restrictive bounds on light neutrino masses, which are however strongly model-dependent. 
We do not discuss this subject here, covered in other contribution to this Special Issue [1]. 

For further details, the reader is referred to recent reviews on neutrino cosmology such as [5l El [7] 
and in particular . A more general review on the connection between particle physics and cosmology 
can be found in [^. 



0.2 The cosmic neutrino background 

The existence of a relic sea of neutrinos is a generic feature of the standard hot big bang model, in num- 
ber only slightly below that of relic photons that constitute the cosmic microwave background (CMB). 
This cosmic neutrino background (CNB) has not been detected yet, but its presence is indirectly es- 
tablished by the accurate agreement between the calculated and observed primordial abundances of 
light elements, as well as from the analysis of the power spectrum of CMB anisotropics and other 
cosmological observables. In this section we will summarize the evolution and main properties of the 
CNB. 



0.2.1 Relic neutrino production and decoupling 

Produced at large temperatures by frequent weak interactions, cosmic neutrinos of any flavour (fe./i.r) 
were kept in equilibrium until these processes became ineffective in the course of the expansion of 
the early Universe. While coupled to the rest of the primeval plasma (relativistic particles such as 
electrons, positrons and photons), neutrinos had a momentum spectrum with an equilibrium Fermi- 
Dirac form with temperature T, 



feqip,T) 



exp 



+ 1 



(1) 



which is just one example of the general case of particles in equilibrium (fermions or bosons, relativistic 
or non- relativistic), as shown e.g. in [lOj. In the previous equation we have included a neutrino chemical 
potential ^i, that would exist in the presence of a neutrino-antineutrino asymmetry, but we will see 



later in section 0.4.1 that even if it exists its contribution can not be very relevant. 

As the Universe cools, the weak interaction rate Ti, falls below the expansion rate and one says 
that neutrinos decouple from the rest of the plasma. An estimate of the decoupling temperature T^ec 
can be found by equating the thermally averaged value of the weak interaction rate 



where cr^ oc Gj^ is the cross section of the electron- neutrino processes with Gp the Fermi constant and 
riu is the neutrino number density, with the expansion rate given by the Hubble parameter H 



Sirp 
3M2 



(3) 
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Here /? oc is the total energy density, dominated by relativistic particles, and Mp = 1 jG^I'^ is the 
Planck mass. If we approximate the numerical factors to unity, with Vy ~ G\T^ and H ~ T'^ jMp^ 
we obtain the rough estimate T^ec ~ 1 MeV. More accurate calculations give slightly higher values of 
Tdec which are flavour dependent because electron neutrinos and antineutrinos are in closer contact 
with electrons and positrons, as shown e.g. in [T]. 

Although neutrino decoupling is not described by a unique Tdec, it can be approximated as an 
instantaneous process. The standard picture of instantaneous neutrino decoupling is very simple 
(see e.g. |10j or and reasonably accurate. In this approximation, the spectrum in eq. ([T]) is 

preserved after decoupling, because both neutrino momenta and temperature redshift identically with 
the expansion of the Universe. In other words, the number density of non-interacting neutrinos remains 
constant in a comoving volume since the decoupling epoch. We will see later that active neutrinos 
cannot possess masses much larger than 1 eV, so they were ultra-relativistic at decoupling. This is 
the reason why the momentum distribution in eq. ([T]) does not depend on the neutrino masses, even 
after decoupling, i.e. there is no neutrino energy in the exponential of feqip)- 

When calculating quantities related to relic neutrinos, one must consider the various possible 
degrees of freedom per flavour. If neutrinos are massless or Majorana particles, there are two degrees 
of freedom for each flavour, one for neutrinos (one negative helicity state) and one for antineutrinos 
(one positive helicity state). Instead, for Dirac neutrinos there are in principle twice more degrees 
of freedom, corresponding to the two helicity states. However, the extra degrees of freedom should 
be included in the computation only if they are populated and brought into equilibrium before the 
time of neutrino decoupling. In practice, the Dirac neutrinos with the "wrong-helicity" states do not 
interact with the plasma at temperatures of the MeV order and have a vanishingly small density with 
respect to the usual left-handed neutrinos (unless neutrinos have masses close to the keV range, as 
explained in section 6.4 of [1], but such a large mass is excluded for active neutrinos). Thus the relic 
density of active neutrinos does not depend on their nature, either Dirac or Majorana particles. 

Shortly after neutrino decoupling the temperature drops below the electron mass, favouring 
annihilations that heat the photons. If one assumes that this entropy transfer did not affect the 
neutrinos because they were already completely decoupled, it is easy to calculate the change in the 
photon temperature before any annihilation and after the electron-positron pairs disappear by 
assuming entropy conservation of the electromagnetic plasma. The result is 



rpaitev / 1 1 \ ^^'^ 



_J_ _ 
J^before ^ ^ 



1.40102, (4) 



which is also the ratio between the temperatures of relic photons and neutrinos T^jTy = (11/4)^/^. 
The evolution of this ratio during the process of annihilations is shown in the left panel of Fig. 
[T| while one can see in the right panel how in this epoch the photon temperature decreases with 
the expansion less than the inverse of the scale factor a. Instead the temperature of the decoupled 
neutrinos always falls as 1/a. 

It turns out that the standard picture of neutrino decoupling described above is slightly modified: 
the processes of neutrino decoupling and annihilations are sufficiently close in time so that some 
relic interactions between and neutrinos exist. These relic processes are more efficient for larger 
neutrino energies, leading to non-thermal distortions in the neutrino spectra at the percent level and 
a slightly smaller increase of the comoving photon temperature, as noted in a series of works listed in 
[U [2] . A proper calculation of the process of non- instantaneous neutrino decoupling demands solving 
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the momentum-dependent Boltzmann equations for the neutrino spectra, a set of integro-differential 
kinetic equations that are difficult to solve numerically. This problem was considered in [12j including 
the effect of flavour neutrino oscillations on the neutrino decoupling process. One finds an increase 
in the neutrino energy densities with respect to the instantaneous decoupling approximation (0.73% 
and 0.52% for z^e's and ly^^s, respectively) and a value of the comoving photon temperature after 
annihilations which is a factor 1.3978 larger, instead of 1.40102. These changes modify the contribution 
of relativistic relic neutrinos to the total energy density which is taken into account using N^s ~ 3.046, 



as defined later in eq. (12). In practice, the distortions calculated in [12j only have small consequences 
on the evolution of cosmological perturbations, and for many purposes they can be safely neglected. 

Any quantity related to relic neutrinos can be calculated after decoupling with the spectrum in 
eq. ([T]) and T^. For instance, the number density per flavour is fixed by the temperature, 

11 ^ " Il7r2 



_ 3 _ 6C(3) 3 

— -I . "-7 — -I . 2 -^7 ' 



which leads to a present value of 113 neutrinos and antineutrinos of each flavour per cm^. Instead, 
the energy density for massive neutrinos should in principle be calculated numerically, with two well- 
defined analytical limits, 

77r2 / 4 \ . 

Puirrii, > Ti,) = myUy . (7) 
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0.2.2 Background evolution 

Let us discuss the evolution of the CNB after decoupling in the expanding Universe, which is described 
by the Friedmann- Robertson- Walker metric 

ds'^ = dt^ - a{tf 6ijdx'dx^ , (8) 

where we assumed negligible spatial curvature. Here a{t) is the scale factor usually normalized to 
unity now (a(to) = 1) and related to the redshift z as a = 1/(1 + z). General relativity tells us the 
relation between the metric and the matter and energy in the Universe via the Einstein equations, 
whose time-time component is the Friedmann equation 

^ = ^^ = ^P = ^o^, (9) 

"J Pc 

that gives the Hubble rate in terms of the total energy density p. At any time, the critical density 
Pc is defined as pc = 3H'^ /SirG, and the current value Hq of the Hubble parameter gives the critical 
density today 

/j° = 1.8788 X 10~2^ h"^ g cm-3 . (10) 
where h = //o/(100 kms^^ Mpc^^). The different contributions to the total energy density are 

P = P-y + Pcdm + Ph + Pu + PA , (11) 

and the evolution of each component is given by the energy conservation law in an expanding Universe 
p = —3H{p + p), where p is the pressure. Thus the homogeneous density of photons p^ scales like 
a~^, that of non-relativistic matter (pcdm for cold dark matter and p^ for baryons) like a~'^, and 
the cosmological constant density p\ is of course time-independent. Instead, the energy density of 
neutrinos contributes to the radiation density at early times but they behave as matter after the 
non-relativistic transition. 

The evolution of all densities is depicted in Fig. [2| starting at MeV temperatures until now. The 
various density fractions = pi/pc are shown in this figure, where it is easy to see which of the 
Universe components is dominant, fixing its expansion rate: first radiation in the form of photons 
and neutrinos (Radiation Domination or RD), then matter which can be CDM, baryons and massive 
neutrinos at late times (Matter Domination or MD) and finally the cosmological constant density 
takes over at low redshift (typically z < 0.5). 

Massive neutrinos are the only particles that present the transition from radiation to matter, when 
their density is clearly enhanced (upper solid lines in Fig. [2]). Obviously the contribution of massive 
neutrinos to the energy density in the non-relativistic limit is a function of the mass (or the sum of 
all masses for which rrii ^ T,^), and the present value fi^ could be of order unity for eV masses (see 
section 0.5 ). 

Shortly after neutrino decoupling, the CNB plays an interesting role in Big Bang Nucleosynthesis 
(BBN), the period of the Universe when the primordial abundances of light elements are created. This 
subject is described in the contribution [3] (for a recent review, see [13] )• Here we just summarize 
the two main effects of relic neutrinos at BBN. The first one is that they contribute to the relativistic 
energy density of the Universe (if rrii, <^ Ti,), thus fixing the expansion rate. This is why BBN gave 
the first allowed range of the number of neutrino species before accelerators (see the next section). 
On the other hand, BBN is the last cosmological epoch sensitive to neutrino flavour, because electron 
neutrinos and antineutrinos play a direct role in the weak processes. We will see some examples of 
BBN bounds on neutrinos (effective number or oscillations) in the following sections. 
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Figure 2: Evolution of the background energy densities in terms of the fractions Qi, from the time 
when Ty = 1 MeV until now, for each component of a flat Universe with h = 0.7 and current density 
fractions JIa = 0.70, fib = 0.05 and ricdm = 1 ~ — ~ ^u- The three neutrino masses are mi = 0, 
m2 = 0.009 eV and jns = 0.05 eV. 



0.3 Extra radiation and the effective number of neutrinos 

Together with photons, in the standard case neutrinos fix the expansion rate during the cosmological 
era when the Universe is dominated by radiation. Their contribution to the total radiation content 
can be parametrized in terms of the effective number of neutrinos A'gff, through the relation 



Pt = + Pv 



7/4x4/3 



P7 > (12) 



where we have normalized to the photon energy density because its value today is known from the 
measurement of the CMB temperature. This equation is valid when neutrino decoupling is complete 
and holds as long as all neutrinos are relativistic. 

We know that the number of light neutrinos sensitive to weak interactions (flavour or active 
neutrinos) equals three from the analysis of the invisible Z-boson width at LEP, Ni, = 2.9840 it 0.0082 
|14j , and we saw in a previous section from the analysis of neutrino decoupling that they contribute as 
iVefj ~ 3.046. Any departure of A'efr from this last value would be due to non-standard neutrino features 
or to the contribution of other relativistic relics. For instance, the energy density of a hypothetical 
scalar particle (j) in equilibrium with the same temperature as neutrinos would be = (7r/30)Tjf, 
leading to a departure of A'gfr from the standard value of 4/7. A detailed discussion of cosmological 
scenarios where A'gfr is not fixed to three can be found in [U [H [T5] . 
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The expansion rate during BBN fixes the produced abundances of hght elements, and in particular 
that of ^He. Thus the value of N^s can be constrained at the BBN epoch from the comparison 
of theoretical predictions and experimental data on the primordial abundances of light elements |31 
I13j . In addition, a value of A'efr different from the standard one would modify the transition epoch 
from a radiation-dominated to a matter-dominated Universe, which has some consequences on some 
cosmological observables such as the power spectrum of CMB anisotropics, leading to independent 
bounds on the radiation content. These are two complementary ways of constraining N^.^ at very 
different epochs. 

A recent analysis of the BBN constraints on N(,q can be found in [16] (see the references therein 
for a list of recent works). The authors have discussed a new and more conservative approach to 
derive BBN constraints on A'efr, motivated by growing concerns on the reliability of astrophysical 
determinations of primordial ^He. According to [16] , the extra radiation at the BBN epoch is limited 
to AA'efT < 1 at 95% C.L. On the other hand, recent analyses of late cosmological observables seem 
to favor N^s > 3, with best-fit values of order 4.3 — 4.4, although with large errorbars as shown 
for instance in |17j . The considered data include CMB temperature anisotropics and polarization, 
combined with other data such as measurements of the present value of the Hubble parameter (Hq), 
the power spectrum of Luminous Red Galaxies (LRG) or distance measurements from the baryon 
acoustic oscillations (BAG) in the distribution of galaxies. The allowed regions at 95% C.L. from [T7] 
are 2.8 < A^eff < 5.9 (WMAP-FBAO-FHo) and 2.7 < N^s < 6.2 (WMAP-hLRG-FHo). Other recent 
analyses of A'eff bounds from cosmological data can be found in [L8^ [19l |20j. These ranges are in 
reasonable agreement with the standard prediction of N^s — 3.046, although they show a marginal 
preference for extra relativistic degrees of freedom whose robustness is still unclear. Moreover, they 
show that there exists an allowed region of N^s values that is common at early (BBN) and more recent 
epochs, although with large errorbars. The upcoming CMB measurements by the Planck satellite 
will soon pin down the radiation content of the Universe, clarifying whether one really needs new 
physics leading to relativistic degrees of freedom beyond the contribution of flavour neutrinos. 



0.4 Neutrino oscillations in the Early Universe 

Nowadays there exist compelling evidences for flavour neutrino oscillations from a variety of experi- 
mental data on solar, atmospheric, reactor and accelerator neutrinos. These are very important results, 
because the existence of flavour change implies that neutrinos mix and have non-zero masses, which 
in turn requires particle physics beyond the Standard Model. Thus it is interesting to check whether 
neutrino oscillations can modify any of the cosmological observables. More on neutrino oscillations 
and their implications can be found in other contributions to this Special Issue or any of the existing 
reviews such as [21 ^ 122 1 l23j. to which we refer the reader for more details. 

It turns out that in the standard cosmological picture all flavour neutrinos were produced with the 



same energy spectrum, as we saw in section 0.2.1 , so we do not expect any effect from the oscillations 
among these three states. This is true up to the small spectral distortion caused by the heating of 
neutrinos from e~^e~ annihilations [12], as described before. In this section we will briefly consider 
two cases where neutrino oscillations could have cosmological consequences: flavour oscillations with 
non-zero relic neutrino asymmetries and active-sterile neutrino oscillations. 
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0.4.1 Active-active neutrino oscillations: relic neutrino asymmetries 

A non-zero relic neutrino asymmetry exists when the number densities of neutrinos and antineutrinos 
of a given flavour are different, and can be parameterized by the ratio r]^^ = {n^^ — njp^)/n^. Such 
a putative asymmetry, that could have been produced by some mechanism well before the thermal 
decoupling epoch, is expected to be of the same order of the cosmological baryon number r]B, i.e. 
a few times 10"^'', from the equilibration of lepton and baryon numbers by sphalerons in the very 
early Universe. In such a case, cosmological neutrino asymmetries would be too small to have any 
observable consequence. However, values fo these parameters which are orders of magnitude larger 
than rjB are not excluded by observations. Actually, large r]u^ are predicted in theoretical models 
where the generation of lepton asymmetry took place after the electroweak phase transition or the 
electroweak washing out of preexisting asymmetries is not effective. 

Neutrino asymmetries can be quantified assuming that a given flavour is characterized by a Fermi- 
Dirac distribution as in eq. with a non-zero chemical potential fi^, or equivalently with the dimen- 
sionless degeneracy parameter = fi^/T (for antineutrinos, = —^u)- In such a case, sometimes 
one says that the relic neutrinos are degenerate (but not in the sense of equal masses). Degenerate 
electron neutrinos have a direct effect on BBN: any change in the z/g/ I'e spectra modifles the primordial 
neutron-to-proton ratio, which in this case is n/p oc exp(— ^j^^). Therefore, a positive decreases the 
primordial ^He mass fraction, while a negative increases it [24j, leading to an allowed range 

- 0.01 < U < 0.07 , (13) 

compatible with = and very restrictive for negative values. In addition a non-zero relic neutrino 
asymmetry always enhances the contribution of the CNB to the relativistic energy density, because 
for any one has a departure from the standard value of the effective number of neutrinos A'efr given, 
if the neutrino spectra follow an equilibrium form, by 



2( 



vr / \ vr 



(14) 



We have seen that this increased radiation modifies the outcome of BBN and that bounds on iVeg can 
be obtained. In addition, another consequence of the extra radiation density is that it postpones the 
epoch of matter-radiation equality, producing observable effects on the spectrum of CMB anisotropies 
and the distribution of cosmic large-scale structures (LSS). Both independent bounds on the radiation 
content can be translated into flavour-independent limits on 

Altogether these cosmological limits on the neutrino chemical potentials or relic neutrino asym- 
metries are not very restrictive, because at least for BBN their effect in the or z^,- sector can be 
compensated by a positive S,^^ . For example, an analysis of the combined effect of a non-zero neutrino 
asymmetry on BBN and CMB/LSS yields the allowed regions [25] 

- O.OK < 0.22, \&,J<2.Q, (15) 

in agreement with similar but more updated bounds as cited in [8] . These limits allow for a very signif- 
icant radiation contribution of degenerate neutrinos, leading many authors to discuss the implications 
of a large neutrino asymmetry in different physical situations (see e.g. [1]). 



It is obvious that the limits in eq. (15) would be modified if neutrino flavour oscillations were 



effective before BBN, changing the initial distribution of flavour asymmetries. Actually, it was shown 
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in |26| [271 [28] that this is the case for the neutrino mixing parameters in the region favoured by 
present data. This resuh is obtained only after the proper inclusion of the refractive terms produced 
by the background neutrinos, which synchronize the oscillations of neutrinos with different momenta 
(which would evolve independently without them). If flavour equilibrium is reached before BBN and 
the momentum distributions of neutrinos keep a Fermi-Dirac form as in eq. Q, the restrictive limits 
on ^j^g in eq. ( |13| ) apply to all flavours. The bounds on the common value of the neutrino degeneracy 
parameter = n^/T would be -0.05 < < 0.07 at 2c7 [29]. 

More recent analyses |30Ll3T] have shown that this conclusion does not always hold. For the present 
measured values of neutrino mixing parameters, the degree of flavour equilibration depends on the 
value of the mixing angle ^13, which fixes the onset of flavour oscillations involving v^s. This in turn 
determines whether neutrinos interact enough with electrons and positrons to transfer the excess of 
energy density due to the initial r]y^ to the electromagnetic plasma. 

In [31] the BBN bounds on the total neutrino asymmetry (cosmological lepton number) were found 
for a range of initial flavour neutrino asymmetries. An example of this analysis is shown in Fig. [3j taken 
from |32) . From this plot one can easily see the effect of flavour oscillations on the BBN constraints 
on the total neutrino asymmetry. With no neutrino mixing the value of r]^^ is severely constrained by 
^He data, while the asymmetry for other neutrino flavours could be much larger. On the other hand, 
flavour oscillations implies that an initially large r/J," can be compensated by an asymmetry in the 
other flavours with opposite sign. The most restrictive BBN bound on r]jj^ applies then to the total 
asymmetry, an effect that can be seen graphically in Fig. [3] as a rotation of the allowed region from 
a quasi-horizontal one for zero mixing to an almost vertical region for sin'^ ^13 = 0.04, in particular 
for the inverted mass hierarchy. In all cases depicted in Fig. [3] the allowed values of the asymmetries 
are fixed by both deuterium and "^He, the latter imposing that the value of rj^^ at BBN must be very 
close to zero. Finally, in [33] it has been shown that '^He data still fix the bounds on rjy when CMB 
results and other cosmological parameters such as neutrino masses are included in the analysis, and 
only for future CMB data the bounds on the asymmetry could be improved. At the same time the 
contribution of neutrino asymmetries to an enhancement of radiation is limited to N^^ < 3.1 — 3.2 for 
values of ^13 in the region allowed by oscillation data. 



0.4.2 Active-sterile neutrino oscillations 

In addition to the flavour or active neutrinos (three species as we saw from accelerator data), there 
could also exist extra massive neutrino states that are sterile, i.e. singlets of the Standard Model gauge 
group and thus insensitive to weak interactions. Most of the current data on neutrino oscillations can 
be perfectly explained with only the three active species, but there exist a few experimental results, 
sometimes called anomalies, that cannot be explained in this framework. If neutrino oscillations are 
responsible for all the experimental data, a solution might require additional, sterile, neutrino species. 
These kind of particles are predicted by many theoretical models beyond the SM, being neutral leptons 
insensitive to weak interactions whose only interaction is gravitational. Their masses are usually heavy, 
while lighter sterile neutrinos are rarer but possible. Here we will briefly summarize these anomalies 
observed in neutrino experiments, whereas an updated review of this subject is given in |34] . 

The long-standing evidence (more than 2>a) for — )• Vf, oscillations comes from the LSND experi- 
ment. Its results pointed out a Am^ ~ 1 eV^, much larger than those required from a three- neutrino 



analysis, as shown later in eq. (21). At the same time, the KARMEN experiment, very similar but 



not identical to LSND, provided no support for such evidence, while a recent data release of the 
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Figure 3: BBN contours at 95% C.L. in the r/i, — r/™ plane for several values of sin^ ^13: (solid line), 
0.04 and normal mass hierarchy (NH) (almost vertical solid line), 0.04 and inverted mass hierarchy 
(IH) (dotted hue) 



MiniBooNE experiment, designed to check the LSND results with larger distance L and energy E 
but similar ratio L/E, is still consistent with the LSND signal. In addition, an unexplained excess 
of electron-like events is observed in MiniBooNE at low energies. The simultaneous interpretation of 
LSND (antineutrino) and MiniBooNE (neutrino and antineutrino) results in terms of sterile neutrino 
oscillations requires CP-violation or some other exotic scenarios, as reviewed in |34) . 

A new anomaly supporting oscillations with sterile neutrinos appeared from a revaluation of reactor 
antineutrino fluxes, which found a 3% increase relative to previous flux calculations. As a result, data 
from reactor neutrino experiments at very short distances can be interpreted as an apparent deficit of 
I'e. This is known as the reactor antineutrino anomaly and is again compatible with sterile neutrinos 
having a Am^ > 1 eV^. Finally, an independent experimental evidence for z/e disappearance at very 
short baselines exists from the Callium radioactive source experiments GALLEX and SAGE. 

The existence of all these experimental hints for sterile neutrinos and a mass scale at the eV is 
intriguing, but so far a fully consistent picture has not emerged. Many analyses have been performed 
trying to explain all experimental data with 1 or 2 additional sterile neutrinos, known as the 3-1-1 or 3-1-2 
schemes, with the corresponding additional mixing parameters. It seems that none of these schemes 
does describe well all data, as explained in detail in [M], but for the topics of this review the potential 
existence of oscillations into sterile neutrinos would lead to important cosmological consequences, such 
as extra radiation from fully or partly thermalized sterile neutrinos or a larger Hot or Warm Dark 
Matter component. The required values of neutrino masses in these 3-1-1 or 3 -|- 2 scenarios are, as 



we will see in section 0.7, in tension with the current cosmological bounds. 

In any case, it is interesting to consider the main effects of additional sterile neutrino species in Cos- 



10 



Neutrino Physics 



Neutrino mass from Cosmology 



mology. Because sterile neutrinos are insensitive to weak interactions, they do not follow the behaviour 
of active neutrinos and are not expected to be present in the early Universe at MeV temperatures. 
Even if they could interact through other kind of interactions, significantly weaker than the standard 
weak ones, as predicted by extended particle physics models, they would have a thermal spectrum 
at very high temperatures but their density would have been strongly diluted by many subsequent 
particle-antiparticle annihilations. Therefore, barring the non-thermal production from additional 
physics beyond the SM, the main way of obtaining a significant abundance of sterile neutrinos is 
through their mixing with the active ones. 

In principle, the cosmological evolution of the active-sterile neutrino system should be found solving 
the corresponding Boltzmann kinetic equations for the density matrices. There exist a vast literature 
on this subject, where different analyses considered several approximations. For a list of references, 
see for instance [H [2]. Although in general one should consider, at least, a 4 x 4 mixing of three 
active neutrinos and one sterile species (with 4 masses, 6 mixing angles and 3 CP-violating phases), 
let us first assume an admixture of one sterile state to electron neutrinos. In the early Universe one 
expects that neutrino oscillations could be effective when the vacuum oscillation term becomes larger 
than the main matter potential term from charged leptons at a temperature Tc- The evolution of the 
active-sterile neutrino system depends on the sign of Am? . For negative values there could be resonant 
oscillations (or resonant production, RP). Instead, for Am^ > one has the so-called non-resonant 
production of sterile neutrinos (NRP). Comparing the value of Tc with the decoupling temperature of 
active neutrinos Tj^d ~ 2 MeV, one obtains that for values 

Am^ < 1.3 X 10^'^eV^ (16) 

active-sterile oscillations are effective after neutrino decoupling. In such a case, the total comoving 
number density of active and sterile neutrinos will be constant because active neutrinos are no longer 
interacting with the rest of the primeval plasma. Correspondingly, distortions in the momentum 
spectra of neutrinos are expected which directly affect the production of ^He at BEN when the active 
neutrinos are of the electron fiavour (see e.g. [35], also for the case of non-zero initial I's abundance). 
Instead, for much larger values of Am^ oscillations are effective before neutrino decoupling, when weak 
interactions are frequent. In this case the actual growth of the sterile neutrino population depends 
on the interplay between oscillations and interactions, while the energy spectrum of the active ones 
will be kept in equilibrium. Their combined contribution to radiation can be as large as N^f^ = 2, 
depending on the specific value of the mixing angle and on the sign of Am^. This "thermalization" of 
the sterile neutrinos is a well-known phenomenon that is very difficult to avoid unless the cosmological 
scenario is drastically modified. 

In the NRP case, it was shown in |36] that the production probability of sterile neutrinos is 

Ts = (sin^ 29m sin^ {^^osct)T a) (17) 

where 9m is the mixing angle in matter and uJqsc the frequency of oscillations in the medium. Here Fq 
is the production rate of active neutrinos in the plasma, and the averaging is made over the thermal 
background. If the oscillation frequency is very high one can substitute sin^(a;osci) — 1/2, and if 9m is 
small one obtains 

Ts « elT, (18) 

for a small number density of Vg and active neutrinos close to equilibrium. Therefore, a thermal or 
close to thermal population of sterile neutrinos is expected provided that such a production rate of 
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Figure 4: Isocountours of the final value of AN^s in the sin^ 29s — ^fn^ plane for vanishing lepton 
asymmetry and > (left panel) and 6m1 < (right pannel). The star denotes the best-fit mixing 
parameters as in the 3+1 global fit in [37]: ((Jm^, sin^ 26's) = (0.9 eV^ 0.089). Adapted from [38] . 
Courtesy of S. Hannestad, I. Tamborra, and T. Tram. 



h's is larger than the expansion rate of the Universe, a condition that holds unless either Am^ or 9 is 
very small. 

For RP of sterile neutrinos at temperatures below Tc the resonance propagates from small to 
large values of the neutrino comoving momentum, covering the whole momentum distribution while 
the active neutrinos are repopulated by interactions. The thermalization of I's is thus significantly 
enhanced, even for quite small values of the mixing angle. 

In order to illustrate this discussion with an actual calculation of the active-sterile system with 
the kinetic equations in the two-flavour approximation, among the many published analyses we have 
chosen a recent one |381. Their results for the final extra contribution of sterile neutrinos to radiation, 
in the case of zero initial lepton asymmetry, are shown as isocontours of AA'eff in Fig. |4] as a function 
of the mixing parameters Jm^ = Ara^ (in eV^ units) and sin^ 29s = sin^ 29. In the NRP case (right 
panel) one can clearly see that the the same final A'cfr corresponds to constant values of sin^29Am?. 
For RP i/s's are more efficiently brought into equilibrium, even for quite small values of the mixing 
angle. All calculations described so far correspond to the two-neutrino limit of one active and one 
sterile states, but a proper calculation should also include the unavoidable presence of mixing among 
the active neutrinos. A full four-flavour calculation has not been performed, but a few analyses did 
solve simplified kinetic equations taking into account active neutrino mixing, such as [39l HOI HI] . 

We note that, for the active-sterile parameters needed to solve the oscillation anomalies described 
at the beginning of this section, the thermalization of sterile neutrinos is achieved, i.e. A'^eff = 4. An 
example of the best-fit values of a particular calculation in the 3 + 1 neutrino model is indicated in 
the plot. Therefore, it seems that such an extra radiation is guaranteed in these situations unless 
oscillations are suppressed, as in the case of a non-zero initial lepton asymmetry r]i/ much larger than 

VB MM- 
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0.5 Massive neutrinos as Dark Matter 

Nowadays the existence of Dark Matter (DM), the dominant non-baryonic component of the matter 
density in the Universe, is well established. A priori, massive neutrinos are excellent DM candidates, in 
particular because we are certain that they exist, in contrast with other candidate particles. Together 
with CMB photons, relic neutrinos can be found anywhere in the Universe with a number density 
given by the present value of eq. ibh of 339 neutrinos and antineutrinos per cm^, and their energy 



density in units of the critical value of the energy density (see eq. (10)) is 



" 93.14 /i2eV ■ ^ 

Here rrii includes all masses of the neutrino states which are non-relativistic today. It is also useful 
to define the neutrino density fraction fu with respect to the total matter density 

= 7 T^-^^ = ^ (20) 

[Pcdm + Ph + Pu) I'm 

In order to check whether relic neutrinos can have a contribution of order unity to the present 
values of fij^ or f^, we should consider which neutrino masses are allowed by non-cosmological data. 
Oscillation experiments measure the differences of squared neutrino masses Amg]^ = — mf and 
Amg]^ = mg — mf, the relevant ones for solar and atmospheric neutrinos, respectively. As a reference, 
we take the following 3a ranges of mixing parameters from }44j (see also [45^ H6] ) , 

Ami,(10-3 eV^) = 2.55t°i^ (-2.43t°:i) 

s?2 = 0.32 ±0.05 

5^3 G [0.36,0.68] ([0.37,0.67]) 

s?3 = 0.0246+|]:|]|J^6 (0.025 ± 0.008) (21) 

Here sfj = sin^ dij, where 0ij {ij = 12, 23 or 13) are the three mixing angles. Unfortunately oscillation 
experiments are insensitive to the absolute scale of neutrino masses, because the knowledge of Am^x > 
and lArn-gj^l leads to the two possible schemes shown in fig. 1 of [8|, but leaves one neutrino mass 
unconstrained. These two schemes are known as normal (NH) and inverted (IH) mass hierarchies, 
characterized by the sign of Am'^i, positive and negative, respectively. In the above equation the 
values in parentheses correspond to the IH, otherwise the mixing parameters present the same allowed 
regions for both hierarchies. For small values of the lightest neutrino mass mo, i.e. mi (ma) for NH 
(IH), the mass states follow a hierarchical scenario, while for masses much larger than the differences 
all neutrinos share in practice the same mass and then we say that they are degenerate. In general, 
the relation between the individual masses and the total neutrino mass can be found numerically, as 
shown in Fig. [5} 

There are two types of laboratory experiments searching for the absolute scale of neutrino masses, a 
crucial piece of information for constructing models of neutrino masses and mixings. The neutrinoless 
double beta decay {Z, A) ^ {Z + 2, A) + 2e~ (in short Oz/2/3) is a rare nuclear processes where lepton 
number is violated and whose observation would mean that neutrinos are Majorana particles. If 
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Figure 5: Allowed values of the total neutrino mass as a function of the lightest state within the 



3(7 regions of the mixing parameters in eq. (21). Blue dotted (red solid) lines correspond to normal 
(inverted) hierarchy for neutrino masses, where mg = mi (mo = Tn^). 



the Oi^2/3 process is mediated by a light neutrino, the results from neutrinoless double beta decay 
experiments are converted into an upper bound or a measurement of the effective mass m^^ 

= |c?2C?3 mi + s?2C?3 ^2 e'^^ + s?3 ^3 e'"^' I , (22) 

where (/>i^2 are the two Majorana phases that appear in lepton-number-violating processes. An impor- 
tant issue for Oz^2/3 results is related to the uncertainties on the corresponding nuclear matrix elements. 
For more details and the current experimental results, see [47| . 

Beta decay experiments, which involve only the kinematics of electrons, are in principle the best 
strategy for measuring directly the neutrino mass [l8]. The current limits from tritium beta decay 
apply only to the range of degenerate neutrino masses, so that ~ niQ, where 

1-2 2 2,22 2,2 2\l/2 /no\ 

= (C12C13 rui + S12C13 m2 + mg) / , (23) 

is the relevant parameter for beta decay experiments. The bound at 95% CL is mo < 2.05 — 2.3 eV 
from the Troitsk and Mainz experiments, respectively. This value is expected to be improved by the 
KATRIN project to reach a discovery potential for 0.3 — 0.35 eV masses (or a sensitivity of 0.2 eV at 
90% CL). Taking into account this upper bound and the minimal values of the total neutrino mass in 
the normal (inverted) hierarchy, the sum of neutrino masses is restricted to the approximate range 

0.06 (0.1) eV < ^ mi < 6 eV (24) 

i 

As we discuss in the next sections, cosmology is at first order sensitive to the total neutrino mass 
S = if states have the same number density, providing information on mo but blind to 
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Figure 6: Allowed regions by oscillation data at the 3a level (eq. (21)) of the three main observables 
sensitive to the absolute scale of neutrino masses. We show the regions in the planes — S and 
mpp — S , see the text for the corresponding definitions. Blue dotted (red solid) regions correspond 
to normal (inverted) hierarchy. 



neutrino mixing angles or possible CP violating phases. Thus cosmological results are complementary 
to terrestrial experiments. The interested reader can find the allowed regions in the parameter space 
defined by any pair of parameters (X]) ^-/S/S) "^^s) [21^ Il9l I50|. The two cases involving ^ are shown 
in Fig. [6} 

Now we can find the possible present values of Qi, in agreement the approximate bounds of eq. 



(24). Note that even if the three neutrinos are non-degenerate in mass, eq. (19) can be safely applied, 
because we know from neutrino oscillation data that at least two of the neutrino states are non- 
relativistic today, because both (Am^i)i/2 ^ 0.05 eV and (Am^J^/^ ^ o.009 eV are larger than the 
temperature Ti, ~ 1.96 K ~ 1.7 x 10^'* eV. If the third neutrino state is very light and still relativistic. 



its relative contribution to Q^, is negligible and eq. (19) remains an excellent approximation of the 



total density. One finds that is restricted to the approximate range 

0.0013 (0.0022) <n^< 0.13 (25) 

where we already included that h ~ 0.7. This applies only to the standard case of three light active 
neutrinos, while in general a cosmological upper bound on has been used since the 1970s to 
constrain the possible values of neutrino masses. For instance, if we demand that neutrinos should 
not be heavy enough to overdose the Universe {Qu < 1), we obtain an upper bound ^ 45 eV 
(again fixing h = 0.7). Moreover, because from present analyses of cosmological data we know that 
the approximate contribution of matter is — 0.3, the neutrino masses should obey the stronger 
bound X] ~ 15 6V. We see that with this simple argument one obtains a bound which is roughly only 
a factor 2 worse than the bound from tritium beta decay, but of course with the caveats that apply 
to any cosmological analysis. In the three-neutrino case, these bounds should be understood in terms 



15 



Neutrino Physics 



Neutrino mass from Cosmology 



Dark matter particles with a large velocity dispersion such as that of neutrinos are called hot 
dark matter (HDM). The role of neutrinos as HDM particles has been widely discussed since the 
1970s, and the reader can find a historical review in [51]. It was realized in the mid-1980s that HDM 
affects the evolution of cosmological perturbations in a particular way: it erases the density contrasts 
on wavelengths smaller than a mass-dependent free-streaming scale. In a universe dominated by 
HDM, this suppression is in contradiction with various observations. For instance, large objects such 
as superclusters of galaxies form first, while smaller structures like clusters and galaxies form via a 
fragmentation process. This top-down scenario is at odds with the fact that galaxies seem older than 
clusters. 

Given the failure of HDM-dominated scenarios, the attention then turned to cold dark matter 
(CDM) candidates, i.e. particles which were non-relativistic at the epoch when the universe became 
matter-dominated, which provided a better agreement with observations. Still in the mid-1990s it 
appeared that a small mixture of HDM in a universe dominated by CDM fitted better the observational 
data on density fluctuations at small scales than a pure CDM model. However, within the presently 
favoured ACDM model dominated at late times by a cosmological constant (or some form of dark 
energy) there is no need for a significant contribution of HDM. Instead, one can use the available 
cosmological data to find how large the neutrino contribution can be, as we will see later. 

Before concluding this section, we would like to mention the case of a sterile neutrino with a mass 
of the order of a few keV's and a very small mixing with the fiavour neutrinos. Such "heavy" neutrinos 
could be produced by active-sterile oscillations but not fully thermalized, so that they could play the 
role of dark matter and replace the usual CDM component. But due to their large thermal velocity 
(slightly smaller than that of active neutrinos), they would behave as Warm Dark Matter and erase 
small-scale cosmological structures. Their mass can be bounded from below using Lyman-a forest 
data from quasar spectra, and from above using X-ray observations. The viability of this scenario is 
currently under careful examination (see e.g. |52]). 

0.6 Effects of neutrino masses on cosmology 

In this section we will briefiy describe the main cosmological observables and the effects that neutrino 
masses cause on them. A more detailed discussion of the effects of massive neutrinos on the evolution 
of cosmological perturbations can be found in Sees. 4.5 and 4.6 of [8]. 

0.6.1 Brief description of cosmological observables 

Although there exist many different types of cosmological measurements, here we will restrict the 
discussion to those that are at present the more important for obtaining an upper bound or eventually 
a measurement of neutrino masses. 

First of all, we have the CMB temperature anisotropy power spectrum, defined as the angular 
two-point correlation function of CMB maps 6T/T{n) (n being a direction in the sky). This function 
is usually expanded in Legendre multipoles 




(26) 
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where Pi{x) are the Legendre polynomials. For Gaussian fluctuations, all the information is encoded 
in the multipoles C/ which probe correlations on angular scales 6 = ir/l. We have seen that each 
neutrino state can only have a mass of the order of 1 eV, so that the transition of relic neutrinos to the 
non-relativistic regime is expected to take place after the time of recombination between electrons and 
nucleons, i.e. after photon decoupling. Because the shape of the CMB spectrum is related mainly to the 
physical evolution before recombination, it will only be marginally affected by the neutrino mass, except 
through a modified background evolution and some secondary anisotropy corrections. There exists 
interesting complementary information to the temperature power spectrum if the CMB polarization 
is measured, and currently we have some less precise data on the temperature x E-polarization (TE) 
correlation function and the E-polarization self-correlation spectrum (EE). 

The current Large Scale Structure (LSS) of the Universe is probed by the matter power spectrum, 
observed with various techniques described in the next section (directly or indirectly, today or in the 
near past at redshift z). It is defined as the two-point correlation function of non-relativistic matter 
fluctuations in Fourier space 

P{k,z) = {\S^{k,z)\') , (27) 

where 6m = Spra/pm- Usually P{k) refers to the matter power spectrum evaluated today (at z = 0). 
In the case of several fluids (e.g. CDM, baryons and non-relativistic neutrinos), the total matter 
perturbation can be expanded as 

Sm = . (28) 

l^t pi 

Because the energy density is related to the mass density of non-relativistic matter through E = mc? , 
(5m represents indifferently the energy or mass power spectrum. The shape of the matter power 
spectrum is affected in a scale-dependent way by the free-streaming caused by small neutrino masses 
of 0(eV) and thus it is the key observable for constraining my with cosmological methods. 



0.6.2 Neutrino free-streaming 

After thermal decoupling, relic neutrinos constitute a collisionless fluid, where the individual particles 
free-stream with a characteristic velocity that, in average, is the thermal velocity t^th- It is possible 
to define a horizon as the typical distance on which particles travel between time ti and t. When the 
Universe was dominated by radiation or matter t ^ ti, this horizon is, as usual, asymptotically equal 
to Vth/H, up to a numerical factor of order one. Similar to the definition of the Jeans length (see 
section 4.4 in [8j), we can define the neutrino free-streaming wavenumber and length as 

As long as neutrinos are relativistic, they travel at the speed of light and their free-streaming length is 
simply equal to the Hubble radius. When they become non-relativistic, their thermal velocity decays 
like 

vth = — - = — - 158(1 + z) kms , (31) 

m m m \ a J \ m I 
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where we used for the present neutrino temperature Tj^ ~ (4/ll)^/^T^ and ~ 2.726 K. This gives 
for the free-streaming wavelength and wavenumber during matter or A domination 



kpsit) = o.8 ''-^J':"'\\^^^ {^r)hMpc-\ (32) 



(l + z)2 VleV 
Xpsit) = 8 , (—] /I'^Mpc , (33) 

where Q.\ and f^m are the cosmological constant and matter density fractions, respectively, evaluated 
today. After the non-relativistic transition and during matter domination, the free-streaming length 
continues to increase, but only like {aH)~^ cx t^^^, i.e. more slowly than the scale factor a cx t^/^. 
Therefore, the comoving free-streaming length Xps/a actually decreases like {a?H)~^ oc t~^^^. As a 
consequence, for neutrinos becoming non-relativistic during matter domination, the comoving free- 
streaming wavenumber passes through a minimum fcnr at the time of the transition, i.e. when m = 
{p) = 3.15Tiy and ao/a = (1 + z) = 2.0 x 10^(m/leV). This minimum value is found to be 

1 /2 

fcnr ^ 0.018 (^^) /iMpc-1 . (34) 

The physical effect of free-streaming is to damp small-scale neutrino density fluctuations: neutrinos 
cannot be confined into (or kept outside of) regions smaller than the free-streaming length, because 
their velocity is greater than the escape velocity from gravitational potential wells on those scales. 
Instead, on scales much larger than the free-streaming scale, the neutrino velocity can be effectively 
considered as vanishing, and after the non-relativistic transition the neutrino perturbations behave 
like CDM perturbations. In particular, modes with k < fcnr are never affected by free-streaming and 
evolve like being in a pure ACDM model. 



0.6.3 Impact of massive neutrinos on the matter power spectrum 

The small initial cosmological perturbations evolve within the linear regime at early times. During 
matter domination, the smallest cosmological scales start evolving non-linearily, leading to the for- 
mation of the structures we see today. In the recent universe, the largest observable scales are still 
related to the linear evolution, while other scales can only be understood using non-linear N-body 
simulations. We will not review here all the details of this complicated evolution (see [H 13 E] and 
references therein) , but we will emphasize the main effects caused by massive neutrinos on linear scales 
in the framework of the standard cosmological scenario: a A Mixed Dark Matter (AMDM) model, 
where Mixed refers to the inclusion of some HDM component. 

On large scales (i.e. on wave-numbers smaller than the value fcnr defined in the previous subsection), 
neutrino free-streaming can be ignored, and neutrino perturbations are indistinguishable from CDM 
perturbations. On those scales, the matter power spectrum P{k, z) can be shown to depend only 
on the matter density fraction today (including neutrinos), and on the primordial perturbation 
spectrum. If the neutrino mass is varied with Q,^ fixed, the large -scale power spectrum remains 
invariant. 

On small scales such that k > k^^, the matter power spectrum is affected by neutrino masses for 
essentially three reasons: 
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1. massive neutrinos do not cluster on those scales. The matter power spectrum can be expanded 
as a function of the three non-relativistic species, 

= ^m^ (l^cdm'^cdm + ^^b<^b + ^^i^^i^l^) . (35) 

On scales of interest and in the recent universe, baryon and CDM fluctuations are equal to 
each other, while 5^ <C (^cdm- Hence, even if the evolution of ^cdm was not affected by neutrino 
masses (which is not the case, see the remaining two points below), the power spectrum would 
be reduced by a factor (1 — fn)"^ with 

U^^- (36) 

2. the redshift of radiation-to-matter equality Zeq or the baryon-to-CDM ratio Wb/^^^cdm might be 
slightly affected by neutrino masses, with a potential impact on the small-scale matter power 
spectrum. This depends very much on which other parameters are kept fixed when the neutrino 
mass varies. If neutrino masses are smaller than roughly 0.5 eV, they are still relativistic at 
the time of radiation-to-matter equality, and the redshift of equality depends on {uj\^ + Wcdm)) 
not on Wm- It is possible to increase and with fixed parameters Q.ra-, ^h-, (^cdm (provided 
that the Hubble parameter also increases like the square root of uj^). In that case, there is 
no significant impact of massive neutrinos on the matter power spectrum through background 
effects, i.e. through a change in the homogeneous cosmological evolution. However, there are 
some important perturbation effects that we will now summarize. 

3. the growth rate of cold dark matter perturbations is reduced through an absence of gravitational 
back-reaction effects from free-streaming neutrinos. This growth rate is set by an equation of 
the type 

'5cdm + ^5cdm = -A;V, (37) 

where (5cdm stands for the CDM relative density perturbation in Fourier space, and ip for the 
metric perturbation playing the role of the Newtonian potential inside the Hubble radius. There 
is a similar equation for decoupled baryons, and very soon after baryon decoupling we can identify 
= <^cdm on scales of interest. On the right-hand side, we neglected time derivatives of metric 
fluctuations playing only a minor role. The right-hand side represents gravitational clustering, 
and is given by the Poisson equation as a function of the total density fluctuation. The second 
term on the left-hand side represents Hubble friction, i.e. the fact that the cosmological expansion 
enhances distances, reduces gravitational forces and slows down gravitational clustering. The 
coefficient a' / a is given by the Priedmann equation as a function of the total background density. 
In a universe such that all species present in the Priedmann equation do cluster, as it is the case 
in a matter-dominated universe with (Jptotal — ^Pcdm + <^Pfe a-nd ptotal = Pcdm + Ph-, the solution is 
simply given by 5cdm oc a: the so-called linear growth factor is proportional to the scale factor. 
But whenever one of the species contributing to the background expansion does not cluster 
efficiently, it can be neglected in the Poisson equation. In that case, the term on the right-hand 
side becomes smaller with respect to the Hubble friction term, and CDM (as well as baryons) 
clusters at a slower rate. This is the case in presence of massive neutrinos and for k ^ k^r- the 
linear growth rate during matter domination is then equal to a^~^^^-^''. During A domination. 



P{k,z) 



SPcdm + Spb + dpp 
Pcdm + Ph+ Pu 
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Figure 7: Ratio of the matter power spectrum including three degenerate massive neutrinos with 
density fraction fu to that with three massless neutrinos. The parameters {uj^, ^}\) = (0.147,0.70) 
are kept fixed, and from top to bottom the curves correspond to f^, = 0.01,0.02,0.03, . . . ,0.10. The 
individual masses rrii, range from 0.046 to 0.46 eV, and the scale knr from 2.1 x 10~^/iMpc^^ to 
6.7 X 10^'^/iMpc^^ as shown on the top of the figure. From [8]. 



this effect sums up with that of the cosmological constant (or of any non-clustering dark energy) , 
which also tends to reduce the growth rate for the very same reason. 

In summary, the small-scale matter power spectrum P{k > fcnr) is reduced in presence of massive 
neutrinos for at least two reason: by the absence of neutrino perturbations in the total matter power 
spectrum, and by a slower growth rate of CDM/baryon perturbations at late times. The third effect 
has the largest amplitude. At low redhsift z ~ 0, the step-like suppression of P{k) starts at A: > k^r and 
saturates at /c ~ 1/i/Mpc with a constant amplitude A.P{k) / P{k) ~ —Sfy. This result was obtained 
by fitting numerical simulations [53], but a more accurate approximation can be derived analytically 
[H [H]. As mentioned in the second item above, neutrino masses can have additional indirect effects 
through a change in the background evolution, depending on which cosmological parameters are kept 
fixed when My varies. 

When fitting data, one can use analytical approximations to the full MDM or AMDM matter power 
spectrum, valid for arbitrary scales and redshifts, as listed in [8|. However, nowadays the analyses are 
performed using the matter power spectra calculated by Boltzmann codes such as CAMB [53j or CLASS 
[55], that solve numerically the evolution of the cosmological perturbations. The step-like suppression 
of the matter power spectrum induced by various values of fu is shown in Fig. [7| 

Is it possible to mimic the effect of massive neutrinos on the matter power spectrum with some 
combination of other cosmological parameters? If so, one would say that a parameter degeneracy exists, 
reducing the sensitivity to neutrino masses. This possibility depends on the intervals [/cmin, fcmax] and 
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[zmin, 2max] in whicli P{k, z) Can be accurately measured. Ideally, if we could have A^min < lO^^/i Mpc^-*^ 
and A;max ^ 1 /iMpc^"*^, the effect of the neutrino mass would be non-degenerate, because of its very 
characteristic step-like effect. Moreover, because neutrinos render the linear growth factor scale- 
dependent, with 5m{k,z) proportional to a (resp. or a}~'^/^f'') for k < k^^ (resp. k ^> knr), the 
amplitude of the step-like suppression is redshift dependent. Using e.g. weak lensing techniques or 
Lyman-a forest data, one could get accurate measurements of the matter spectrum at at several 
redshifts in the range < 2; < 3. This will offer an opportunity to test the "time dependence" of the 
neutrino mass effect, and to distinguish it, for instance, from an equivalent step-like suppression in 
the primordial spectrum, that would still imply a scale-independent growth factor. 

The problem is that usually the matter power spectrum can only be accurately measured in 
the intermediate region where the mass effect is neither null nor maximal: in other words, many 
experiments only have access to the transition region in the step-like transfer function. In this region, 
the neutrino mass affects the slope of the matter power spectrum in a way which can be easily 
confused with the effect of other cosmological parameters. Because of these parameter degeneracies, 
the current LSS data alone cannot provide significant constraints on the neutrino mass, and it is 
necessary to combine them with other cosmological data, in particular the CMB anisotropy spectrum, 
which could lift most of the degeneracies. Still, for exotic models with e.g. extra relativistic degrees of 
freedom, a constant equation-of-state parameter of the dark energy different from —1 or a non-power- 
law primordial spectrum, the neutrino mass bound can become significantly weaker. 

The value of fcmax is not limited by observational sensitivities, but by the range in which we trust 
predictions from linear theory. Beyond the scale of non-linearity (of the order of fcmax = 0.15/iMpc~^ 
at 2; = 0), the data is only useful provided that one is able to make accurate predictions not only for 
the non-linear power spectrum, but also for redshift-space distortions (coming from the fact that we 
observe the redshift of objects, not their distance away from us), and finally for the scale dependence of 
the light-to- mass bias (relating the power spectrum of a given category of observed compact objects to 
the underlying total matter power spectrum). Spectacular progresses are being carried on these three 
fronts, thanks to better N-body simulations and analytical techniques. Including massive neutrinos in 
such calculations is particularly difficult, but successful attempts were presented e.g. in [56l \57\ l58]. 
The neutrino mass impact on the non-linear matter power spectrum is now modeled with rather good 
precision, at least within one decade above A;max in wave-number space. It appears that the step-like 
suppression is enhanced by non-linear effects up to roughly AP{k)/P{k) ~ —10/^ (at redshift zero 
and /c ~ 1 /iMpc~^), and is reduced above this scale. Hence, non-linear corrections render the neutrino 
mass effect even more characteristic than in the linear theory, and may help to increase the sensitivity 
of future experiments. 

Until this point, we reduced the neutrino mass effect to that of the parameter fiy or M^. In 
principle, the mass splitting between the three different families for a common total mass is visible 
in the matter power spectrum. The time at which each species becomes non-relativistic depends on 
individual masses mj. Hence, both the scale of the step-like suppression induced by each neutrino and 
the amount of suppression in the small-scale power spectrum have a small dependence on individual 
masses. The differences between the power spectrum of various models with the same total mass and 
different mass splittings was computed numerically in [59j for the linear spectrum, and [58j for the 
non-linear spectrum. At the moment, it seems that even the most ambitious future surveys will not 
be able to distinguish these mass splitting effects with a reasonable significance |60^ 161] . 
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0.6.4 Impact of massive neutrinos on the CMB anisotropy spectrum 

For neutrino masses of the order of 1 eV (about fi, < 0.1) the three neutrino species are still relativistic 
at the time of photon decoupling, and the direct effect of free-streaming neutrinos on the evolution of 
the baryon-photon acoustic oscillations is the same in the ACDM and AMDM cases. Therefore, the 
effect of the mass can only appear at two levels: that of the background evolution, and that of secondary 
CMB anisotropics, related to the behavior of photon perturbations after decoupling. Both levels are 
potentially affected by the evolution of neutrinos after the time of their non-relativistic transition. If 
neutrinos were heavier than a few eV, they would already be non-relativistic at decoupling, and they 
could trigger more direct effects in the CMB, as described in [62]. However, we will see later that this 
situation is disfavoured by current upper bounds on the neutrino mass. 

Let us first review the background effects of massive neutrinos on the CMB. Because the temper- 
ature and polarization spectrum shape is the result of several intricate effects, one cannot discuss the 
neutrino mass impact without specifying which other parameters are kept fixed. Neutrinos with a 
mass in the range from 10~^ eV to 1 eV should be counted as radiation at the time of equality, and as 
non-relativistic matter today: the total non-relativistic density, parametrized by uj^n = ^mh"^, depends 
on the total neutrino mass M^, = ^ ■ mj. Hence, when My is varied, there must be a variation either 
in the redshift of matter-to-radiation equality Zeq, or in the matter density today ojm- 

This can potentially impact the CMB in three ways. A shift in the redshift of equality affects the 
position and amplitude of the peaks. A change in the non-relativistic matter density at late times 
can impact both the angular diameter distance to the last scattering surface dA{zdcc), controlling 
the overall position of CMB spectrum features in multipole space, and the slope of the low-/ tail of 
the CMB spectrum, due to the late Integrated Sachs- Wolfe (ISW) effect. Out of these three effects 
(changes in Zeq, in iIa and in the late ISW), only two can be cancelled by a simultaneous variation of 
the total neutrino mass and of other free parameters in the AMDM model. Hence, the CMB spectrum 
is sensitive to the background effect of the total neutrino mass. In practice however, the late ISW effect 
is difficult to measure due to cosmic variance and CMB data alone cannot provide a useful information 
on sub-eV neutrino masses. If one considers extensions of the AMDM, this becomes even more clear: 
by playing with the spatial curvature, one can neutralize all three effects simultaneously. But as soon 
as CMB data is used in combination with other background cosmology observations (constraining for 
instance the Hubble parameter, the cosmological constant value or the BAO scale), some bounds can 
be derived on My. 

There exists another effect of massive neutrinos on the CMB at the level of secondary anisotropics: 
when neutrinos become non-relativistic, they reduce the time variation of the gravitational potential 
inside the Hubble radius. This affects the photon temperature through the early ISW effect, and 
leads to a depletion in the temperature spectrum of the order of (ACi/Ci) ~ —{M^/0.1eV)% on 
multipoles 20 < / < 500, with a dependence of the maximum depletion scale on individual masses mj. 
This effect is roughly ten times smaller than the depletion in the small-scale matter power spectrum, 
AP{k)/P{k) ~ -(M^/0.01eV)%. 

We show in Figure [8] the effect on the CMB temperature spectrum of increasing the neutrino mass 
while keeping Zeq and fixed: the only observed differences are then for 2 < / < 50 (late ISW 
effect due to neutrino background evolution) and for 50 < / < 200 (early ISW effect due to neutrino 
perturbations). 

We conclude that the CMB alone is not a very powerful tool for constraining sub-eV neutrino 
masses, and should be used in combination with homogeneous cosmology constraints and/or measure- 
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Figure 8: CMB temperature spectrum with different neutrino masses. Some of the parameters of the 
AMDM model have been varied together with M,y in order to keep fixed the redshift of equality and 
the angular diameter distance to last scattering. 



ments of the LSS power spectrum, for instance from galaxy clustering, galaxy lensing or CMB lensing 



(see section 0.8). 



0.7 Current bounds on neutrino masses 

In this section we review how the available cosmological data are used to get information on the 
absolute scale of neutrino masses, complementary to laboratory experiments. Note that the bounds in 
the next subsections are all based on the Bayesian inference method, and the upper bounds on the sum 
of neutrino masses are given at 95% C.L. after marginalization over all free cosmological parameters. 
We refer the reader to section 5.1 of |8J for a detailed discussion on this statistical method, as well as 
for most of the references for the experimental data or parameter analysis. 

0.7.1 CMB anisotropy bounds 

The experimental situation of the measurement of the CMB anisotropics is dominated by the seven- 
year release of WMAP data pJT], which improved the already precise TT and TE angular power spectra 
of the previous releases, and included a detection of the E-polarization self-correlation spectrum (EE). 
On similar or smaller angular scales than WMAP, we have results from experiments that are either 
ground-based or balloon-borne (ACT, SPT, . . . ). 

We saw in the previous section that the CMB spectrum has a small sensitivity to neutrino masses 
even when each mass is below 0.5 eV, and all non-relativistic transitions take place after photon 
decoupling [HS]. This sensitivity is due to background effects (mainly the late ISW effect if all other 
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AMDM parameters are left free) plus perturbation effects (mainly the early ISW effect). Therefore, it 
is possible to constrain neutrino masses using CMB experiments only, down to the level at which these 
small effects are masked by instrumental noise, or by cosmic variance, or by parameter degeneracies in 
the case of some cosmological models beyond the minimal A Mixed Dark matter framework. WMAP 
alone is able to set a limit < 1.3 — 1.4 eV depending on whether the dark energy component is a 
cosmological constant or not [T7j. Because the neutrino mass effects in the CMB are visible mainly at 
I < 500, combining WMAP data with other CMB datasets (from e.g. ACT or SPT) does not improve 
this bound. On the contrary, adding more information on the background cosmological evolution is 
very useful, because it helps removing degeneracies between the various parameters: in that case, 
the CMB data can better probe neutrino masses through their background effects. For instance, the 
bound from WMAP combined with BAO scale measurements and a direct determination of Hq by 
the Hubble space telescope Key Project [64j is significantly stronger: < 0.58 eV (95% CL) [17], 
while the combination of WMAP with a different Hq determination at various redshift from early-type 
galaxy evolution gives Mi, < 0.48 eV (95% CL) [53]. These are important results, because they do 
not depend on the uncertainties from LSS data discussed next. 

0.7.2 Large Scale Structure observations 

The matter power spectrum can be probed with various methods on different scales and redshifts. Let 
us review here the major techniques which have led so far to relevant neutrino mass bounds. 

Galaxy power spectrum. Galaxy maps (or, similarly, cluster maps) can be smoothed over 
small scales and Fourier transformed in order to provide a power spectrum. Relating such a spectrum 
to the total matter power spectrum is a tricky exercise, especially on small scales (corresponding to 
wave-numbers k > 0.1/iMpc~^), because it is difficult to make accurate predictions for non-linear 
corrections, redshift-space distortions and the light-to-mass bias. Current bounds are based on the 
analysis of linear scales only, for which the assumption of a scale-independent bias is well motivated. 
This bias is however left as a free parameter, in such way that galaxy spectrum data give indications on 
the shape, but not on the global amplitude of P{k, z). In the next subsection, we will report constraints 
from the halo power spectrum of Large Red Galaxies (denoted later as Gal-LRG), measured by [66] . 
the spectrum of the MegaZ catalogue (Gal-MegaZ), used by [HT], and the WiggleZ Dark Energy Survey 
(Gal-WiggleZ), used by [68]. The first two data sets are actually extracted from the same big survey, 
the Sloan Digital Sky Survey (SDSS). 

For sufficiently deep galaxy surveys, it is possible to separate galaxies into redshift bins and compute 
different correlation functions at different redshifts. This technique, called tomography, can be very 
useful for constraining the scale-dependent growth factor induced by neutrino masses. In that case, 
the data can be reduced to a set of two-dimension power spectra in different shells, each of them 
related to P{k, z) in a narrow redshift range. Recently, such a tomographic analysis was used by [69] 
for constraining neutrinos, using galaxies from the Canada-France-Hawaii- Telescope Legacy Survey 
(CFHTLS), split in three redshift bins covering the ranges 0.5 < z < 0.6, 0.6 < z < 0.8 and 0.8 < z < 
1.0 (this dataset will be denoted as Gal-CFHTLS). 

Cluster mass function. Instead of probing directly the matter power spectrum P{k, z) from the 
spatial distribution of objects, it is possible to constrain integrated quantities of the type / dk P{k)W {k), 
where W{k) stands for a given window function. One such quantity is related to the histogram of 
cluster masses. If the mass of a significant number of galaxy clusters within a given redshift bin is 
known, this histogram gives an estimate of the so-called cluster mass function, dn(M , z) / dM , with 
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dn being the number of clusters of redshift approximately equal to z, and with a mass in the range 
[M, M + dM]. This function is related to a^{M, z), the variance of the density in spheres enclosing a 
mass M, itself derived from the convolution of the power spectrum P{k, z) with an appropriate window 
function. In the next subsection we will refer to bounds derived from cluster abundances probed by 
X-ray observations from the ROSAT survey, presented by [70] (denoted later as Clus-ROSAT), and by 
optical observations from the MaxBCG catalogue, presented by [H] (denoted later as Clus-MaxBCG). 

Galaxy weak lensing. The image of observed galaxies is distorted by gravitational lensing 
effects, caused by density fluctuations along the line of sight. One of these effects is called cosmic 
shear. It corresponds to the squeezing of an image in one direction in the sky, and its stretching in 
the orthogonal direction. Because such distortions are coherent over the angular size of the lensing 
potential wells responsible for lensing, they tend to align slightly the apparent major axis of galaxies 
in a given patch of the sky. Hence the average cosmic shear in a given direction can be estimated 
statistically by averaging over major axis orientations. The analysis of a catalogue of images leads 
to a map of the lensing potential, itself related to the matter power spectrum P{k, z) through the 
Poisson equation. If the number of source galaxies is sufficient, it is possible to split the catalogue in 
several redshift bins, and to obtain a three-dimensional reconstruction of the gravitational potential 
in our past-line cone, and of P{k, z) at various redshifts. Cosmic shear tomography is particularly 
useful for measuring neutrino masses, because it can probe the scale-dependent growth factor induced 
by neutrino masses over an extended range of redshifts. Current cosmic shear surveys already allow 
to put bounds on neutrino parameters: we will refer later to results from the CFHTLS presented by 
[72] and denoted as WL-CFHTLS. 

Lyman alpha forests. The most luminous and distant compact objects that we can observe are 
quasars. Some of the photons emitted by quasars interact along the line-of-sight. In particular, a 
fraction of photons are absorbed at the Lyman alpha wavelength by hydrogen atoms located in the 
Interstellar Galactic Medium (IGM). The absorbed fraction in a given point of the photon trajectory 
is proportional to the local density of neutral hydrogen. Because photons are continuously redshifted, 
absorption in a given point is seen by the observer as a depletion of the spectrum at a given frequency. 
Hence, inside a limited range called the Lyman alpha forest, the frequency-dependence of quasar 
spectra is a tracer of the spatial fluctuations of the hydrogen density along the line of sight. Lyman 
alpha forests in quasar spectra offer an opportunity to reconstruct the hydrogen density fluctuation 
along several line-of-sights in a given redshift range. After Fourier expanding each spectrum and 
averaging over many spectra, one gets an estimate of the so-called flux power spectrum Pf{z, k), that 
can be related to the total matter power spectrum P{k,z). Unfortunately, the flux power spectrum 
does not probe linear scales, but mildly non-linear scales. In order to relate Pp(z,A;) to P{k), it is 
necessary to perform N-body simulations with a hydrodynamical treatment of baryons, accounting 
for the complicated thermodynamical evolution of the IGM (which depends on star formation). Also, 
a limitation of this technique comes from the fact that the emitted quasar spectra already have a 
non-trivial frequency dependence, and that photons are affected by several other effects than Lyman 
alpha absorption along the line-of-sight. Nevertheless, a careful modeling of all relevant effects allows 
to obtain interesting constraints. The fact that Lyman alpha forests probe mildly non-linear scales 
rather than strongly non-linear ones is of course crucial for keeping systematic errors under control. 
Lyman alpha observations typically constrain the matter power spectrum in the wavenumber range 
0.3 < k < 3 h/Mpc and in the redshift range 2 < z < 5. We mention below some neutrino mass bounds 
inferred from quasar spectra obtained by the SDSS and presented in [73], denoted as Ly-a-SDSS. 
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Cosmological data 


Reference 


w = —1 


w 7= — i 


W MAr^ ( -+-Kji8iL-Lit\\j-\-iiQ 




n /I /I 
U.44 


U. / 


WMAP5+Gal-MegaZ+BA0+SNIa 


[67J 


0.325 


0.491 


WMAP5+Gal-MegaZ+BAO+SNIa+i?o 


167J 


0.281 


0.471 


WMAP7+Gal-WiggleZ 


[68J 
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0.29 
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0.64(0.44) 
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[69J 


0.41(0.29) 
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[70] 


0.33 


0.43 


WMAP5+ifo+Clus-MaxBCG 


[71J 


0.40 


0.47 


WMAP5+BA0+SNIa+WL-CFHTLS 


[72J 


0.53 





Table 1: 95% CL upper bounds on the total neutrino mass M^, in eV, for various combinations of 
CMB, homogeneous cosmology and LSS data sets. The first seven lines refer to galaxy power spectrum 
measurements, the next two lines to cluster mass function measurements, the last line to a weak lensing 
survey. WMAP5 or 7 stands for WMAP 5 or 7- year data. Hq refers to the direct measurement by the 
HST Key Project [M], and BAO to estimates of the scale of Baryon Acoustic Oscillations at various 
redshifts. Other acronyms refer to various Large Scale Structure dataset referred in section 0.7.2, In 
the last column, the cosmological constant was replaced by a Dark Energy component with arbitrary 
equation-of-state parameter w. 



0.7.3 Large Scale Structure bounds 

Using LSS observations in combination with CMB data offers an opportunity to observe (or to bound) 
the step-like suppression of the matter power spectrum in presence of neutrino masses, as explained in 
section 0.6.3 and illustrated in Fig. [7j The use of CMB data is crucial in order to constrain parameters 



like the baryon density, the primordial spectrum amplitude, the tilt, and a combination of ujm and h. 
Without such constraints, there would be too much freedom in the matter power spectrum fitted to 
LSS data for identifying a smooth step-like suppression. 

We summarize in Table [T] the main constraints on Mi, available at the time of writing, obtained 
from combinations of CMB plus homogeneous cosmology data, galaxy power spectrum data and 
cluster abundance data. Current data set are far from reaching the sensitivity required to probe the 
mass splitting of the total mass M^y = niu- between different species. The constraints mentioned 
below have been derived in the case of three degenerate neutrinos with mass rrii, = but they 

roughly apply to the total mass of any scenario. Also, the bounds shown in the first column of 
Table [T] have been obtained under the assumption of a minimal ACDM model with massive neutrinos, 
featuring seven free parameters. More conservative bounds are sometimes derived for basic extensions 
of this model, with one or two more parameters. The constraints do not change significantly when 
assuming, for instance, a primordial spectrum with a running of the tilt [din Ug/ dink], or a significant 
contribution to the CMB of primordial gravitational waves [7l]. Parameters known to be slightly 
degenerate with neutrino masses and leading to weaker bounds are w, the equation-of-state parameter 
of a Dark Energy component (substituting the cosmological constant), and A^eff; the effective number 



of neutrinos discussed in section 0.3, The degeneracy with w, explained in [74], is illustrated by the 



last column in Table [T] We did not include in this table current limits from Lyman alpha forest 
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data: this is a dehcate matter for the reasons mentioned previously, and a careful modeling of all 
systematic effects leads to rather weak neutrino mass bounds. The conservative analysis of |73j, based 
on Ly-a-SDSS, gives a bound M,j < 0.9 eV (95% CL) from Ly-a-SDSS data alone. 

In conclusion, the combination of available data sets consistently indicates that the total neutrino 
mass is below 0.3 eV at the 95% CL (0.5 eV if we allow for Dark Energy with arbitrary w). This 
means that the "degenerate scenario" in which all neutrinos share roughly the same mass is almost 
excluded. The data is about to probe the region in which masses are different from each other, and 
are ranked according to the NH or IH scenario. Other recent summaries of existing bounds have been 
recently presented in P [T71 [5011711 175] . 

0.8 Future sensitivities on neutrino masses from cosmology 

Future CMB observations will have increasing sensitivity to neutrino masses, not only thanks to 
smaller error bars on the temperature and polarization spectra. They will make it possible to probe 
the large scale structure of the universe with a new technique: CMB lensing extraction. The weak 
lensing of the last scattering surface by nearby galaxy clusters induces specific non-gaussian patterns 
in CMB maps that can be extracted using some non-linear estimators. This method allows to measure 
the lensing potential up to z ~ 3, and to infer constraints on P{k, z) at such high redshifts. 

In the very close future, significant improvements on the neutrino mass bounds will be triggered 
by the Planck CMB satellite. The forecasts presented in j76j predict a neutrino mass sensitivity of 
cr(Mi/) ~ 0.1 eV from Planck alone, using the lensing extraction technique of [77J. This would be twice 
better than without lensing extraction. 

Several galaxy surveys with better sensitivity and larger volume are about to release data or have 
been planned over the next decades, including the Baryon Oscillation Spectroscopic Surve}j^ (BOSS), 
the Dark Energy Survej|^(DES), the Large Synoptic Survey Telescop^(LSST) or the Euclid satellit^ 
Also, in [78j it was pointed out that in the future accurate measurements could be inferred from cluster 
surveys. Because clusters are more luminous than galaxies, they can be mapped up to higher redshift. 
Concerning cosmic shear surveys, spectacular improvements are expected from Pan-STARR^ or the 
DES, LSST and the Euclid surveys already mentioned above. 

In a near future, the prediction of ref. |79j is that the combination of Planck (with lensing extrac- 
tion) with BAO scale information from BOSS could lower the error down to a{My) ~ 0.06 eV. In 
addition, the authors of |80] find that adding Lyman alpha data from BOSS should lead to comparable 
sensitivities, and even better results might be expected from the addition of galaxy power spectrum 
data from the same survey. 

With better tomographic data (for either galaxy clustering or cosmic shear), it will become possible 
to probe the scale-dependence of the growth factor induced by neutrino masses (or in other words, the 
fact that the step- like suppression has an amplitude increasing with time), and to reach spectacular 
sensitivites. We present below the typical sensitivities expected for a collection of planned surveys (not 
all approved). These numbers should be taken with care because forecasts are based on an idealization 

^ cosmology . Ibl . go v/BOSS/ 
■^www . darkenergysurvey . org/ 
•^www. Isst . org/ 
'*sci . esa. int/euclid 
^pan-Starrs . if a.hawaii . edu/ 
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of each experiment, as weh as on several assumptions hke the underlying cosmological model, or even 
the fiducial value of the neutrino mass itself. 

In ref. |81] it was found that the measurement of the galaxy harmonic power spectrum in seven 
redshift bins by the DES should lead to a sensitivity of (j{My) ~ 0.06 eV when combined with Planck 
data (without lensing extraction). Similar bounds were derived in |82j for another combination of 
comparable experiments. This shows that at the horizon of 2014 or 2015, a total neutrino mass close 
to Mjj ~ 0.1 eV could be marginally detected at the 2-a level by cosmological observations. Because 
this value coincides with the lowest possible total mass in the inverted hierarchy scenario, the latter 
could start to be marginally ruled out in case the data still prefers My = 0. 

The sensitivity of cosmic shear data from a satellite experiment comparable to Euclid was calculated 
in [83], where it was found that it would shrink to a{My) ~ 0.03 eV in combination with Planck data 
(without lensing extraction). The forecast of [84| based on galaxy clustering data also from Euclid 
(completed at small redhsift by similar data from BOSS) gives comparable numbers. Constraints 
based on the ground-based Large Synoptic Survey Telescope should be slightly weaker [SS]. Hence, in 
the early 2020's, we expect that a combination of cosmological data sets could detect the total neutrino 
mass of the normal hierachy scenario, ~ 0.05 eV, at the 2-a level. If the total mass is instead 
close to My ~ 0.1 eV, it will be detected at the A-a level. However, in that case, available experiments 
would not have enough sensitivity for making the difference between an inverted and normal hierarchy 
scenario with the same M^. 

Even more progress could be provided by the promising technique of 21-cm surveys. Instead of 
mapping the distribution of hydrogen atoms trough the absorption rate of photons traveling from 
quasars, it should be possible to observe directly the photons emitted by these atoms at a wavelength 
A ~ 21 cm from the transition from one hyperfine level to the other. While travelling towards 
the observer, these photons are redshifted, and seen with a wavelength indicating the position of 
the emitting atoms in redshift space. Recent theoretical progresses in this field show that using 
this technique, future dedicated experiments should be able to map hydrogen and hence baryonic 
fluctuations at very high redshift (typically 6 < z < 12), and to probe the matter power spectrum 
deep inside the matter-dominated regime on linear scales |86j . This field is still in its infancy, and the 
forecasts presented so far have to be taken with care, due to the difficulty to make a realistic estimate 
of systematic errors in future data sets. A sensitivity of a{Mu) ~ 0.075 eV for the combination of 
Planck with the Square Kilometer Array (SKA) project, or a{My) ~ 0.0075 eV with the Fast Fourier 
Transform Telescope (FFTT), was found in [6T]. However, the authors show that such impressive 
experiments would still fail in discriminating between the NH and IH scenario. 

An eventual post-Planck CMB satellite or post-Euclid survey would also have a great potential. 
The forecast analysis in [HTj shows that for a CMB satellite of next generation one could get a{My) ~ 
0.03 eV alone, thanks to a very precise reconstruction of the CMB lensing potential, while [78j discusses 
the potential of cluster surveys. Finally, the authors of |60] show how far the characteristics of an 
hypothetical galaxy or cosmic shear survey should be pushed in order to discriminate between two 
allowed NH and IH scenarios with the same total mass. 

0.9 Conclusions 

Neutrinos, despite the weakness of their interactions and their small masses, can play an important role 
in Cosmology that we have reviewed in this contribution. In addition, cosmological data can be used 
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to constrain neutrino properties, providing information on these elusive particles that complements 
the efforts of laboratory experiments. In particular, the data on cosmological observables have been 
used to bound the radiation content of the Universe via the effective number of neutrinos, including 
a potential extra contribution from other relativistic particles. 

But probably the most important contribution of Cosmology to our knowledge of neutrino prop- 
erties is the information it can provide on the absolute scale of neutrino masses. We have seen that 
the analysis of cosmological data can lead to either a bound or a measurement of the sum of neutrino 
masses, an important result complementary to terrestrial experiments such as tritium beta decay and 
neutrinoless double beta decay experiments. In the next future, thanks to the data from new cosmo- 
logical experiments we could even hope to test the minimal values of neutrino masses guaranteed by 
the present evidences for flavour neutrino oscillations. For this and many other reasons, we expect 
that neutrino cosmology will remain an active research field in the next years. 
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